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Abstract.  The  convergence  of  advective  and  diffusive  buoyancy  flux  must 
match  the  air-sea  buoyancy  flux  between  two  outcropping  isopycnals.  This 
leads  to  a diagnostic  framework  for  water  mass  transformation  in  which  a 
myriad  of  different  processes  can  be  incorporated  under  a unifying  balance.  We 
review  how  the  diapycnal  advection  due  to  ubiquitous  mixed  layer  entrainment 
can  be  included  in  this  framework,  and  we  estimate  its  contribution  to  the 
large  scale  transformation.  We  also  consider  how  decomposing  the  flow  and 
buoyancy  field  into  mean,  eddy  and  turbulent  parts  leads  to  clarifying  the 
interaction  of  mixed  layer  and  mesoscale  (or  sub-mesoscale)  eddies  in  the 
overall  large  scale  balance. 


1.  Introduction 

Quantifying  the  water  mass  changes  and  finding  the 
water  mass  formation  rates  for  specific  density  classes 
are  central  to  understanding  the  diabatic  ocean  circula- 
tion. Accurate  quantification  of  water  mass  formation 
rates  is  also  important  for  understanding  variability  in 
oceanic  heat  and  freshwater  fluxes  and  their  contribu- 
tion to  climate. 

The  current  description  of  water  mass  transforma- 
tion was  formulated  by  Walin  (1982),  who  proposed 
that  the  poleward  surface  drift  in  the  ocean  can  be  de- 
termined directly  by  the  air-sea  heat  flux  —Bq  at  the 
surface.1  This  led  Walin  to  suggest  the  relation 

F(T)  = B0 % (1) 

where  F is  the  cross-isotherm  advection  (volumetric 
flow  rate),  and  SS  is  the  area  enclosed  between  the  two 
isotherm  outcrops  as  shown  in  Figure  1. 

This  relationship  was  subsequently  successfully  used 
in  diagnostic  calculations  ( Tziperman  1986,  Speer  and 
Tziperman  1992).  These  calculations  showed  close  cor- 
respondence between  the  water  mass  formation  derived 
by  the  diathermal  advection  from  the  Walin  relation 
and  the  estimated  rates  of  water  mass  formation  for 


1The  interest  in  water  mass  formation  rate  in  specific  density 
classes  has  followed  various  diagnostics,  e.g.  in  Walin  (1982) 
arising  out  of  considering  budgets  with  heat  fluxes  and  isotherms, 
and  in  other  studies,  the  distribution  of  isopycnals  and  buoyancy 
fluxes.  Thus  we  retain  the  symbol  Bo  for  both  surface  heat  fluxes 
and  surface  buoyancy  fluxes. 


Figure  1.  A schematic  of  the  diapycnal  advection  F(T) 
considered  by  Walin  (1982).  The  net  surface  heat  flux 
is  balanced  by  a net  advective  flux  at  the  surface. 


Sub- Tropical  mode  water  and  the  Sub-Polar  Mode  Wa- 
ter from  observations.  However,  the  Walin  (1982)  for- 
mulation assumed  negligible  mixing  in  the  upper  ocean. 

In  the  same  year  that  Walin  proposed  his  ideas  about 
this  diathermal  advection,  Niiler  and  Stevenson  (1982) 
sought  to  constrain  the  values  of  diapycnal  mixing  by 
considering  the  heat  budget  of  the  isotherms  that  out- 
crop in  the  tropics  and  are  subject  to  net  heat  gain 
through  the  year  as  shown  in  Figure  2. 

In  this  case,  since  by  definition  the  net  change  in  the 
mean  annual  isothermal  position  is  zero,  there  cannot 
be  a net  diathermal  advection  and  the  net  heat  flux 
must  be  balanced  by  a net  diffusive  flux  due  to  mix- 
ing. Therefore,  at  least  in  the  limiting  case  of  closed 
mean  isopycnal  surfaces,  the  role  of  diffusive  flux  due 
to  mixing  cannot  be  ignored,  and  the  surface  drift  of 
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Figure  2.  A schematic  of  the  outcropping  tropical  an- 
nual mean  isotherms  considered  by  Niiler  and  Steven- 
son (1982).  The  net  surface  heat  flux  must  be  balanced 
by  a net  diffusive  flux  due  to  mixing. 

Bo 


\ 


Figure  3.  Air-sea  buoyancy  flux  Bo  SS  between  two 
outcropping  isopycnals  b and  b + 5b  must  be  balanced 
by  convergence  of  diapycnal  advective  flux  due  to  ad- 
vection  A(b)  and  diffusive  flux  D(b).  From  Garrett  et 
al.  (1995). 


Walin  (1982)  cannot  achieve  this  balance. 

This  realization  led  Garrett  et  al.  (1995)  to  clarify 
the  Walin  formulation.  They  insisted  that  considering 
the  diapycnal  advection  about  any  averaged  isopycnal 
position  would  require  diapycnal  mixing,  since  other- 
wise the  isopycnal  would  simply  be  advected  around  by 
the  flow.  As  Figure  3 shows,  the  convergence  of  the  di- 
apycnal advective  and  diffusive  fluxes  must  match  the 
air-sea  buoyancy  flux  between  two  isopycnals. 

The  Walin  (1982)  approach  was  thus  generalized  by 


Figure  4.  Garrett  et  al.  (1995)’s  ideal  case.  With  no 
horizontal  mixing  and  only  vertical  mixing,  D(b)  can 
be  ignored  and  the  diapynal  advection  equals  the  air- 
sea  transformation,  or,  A — F. 

Garrett  et  al.  (1995)  as 

A(b)  = F(b)-^  (2) 

where  D(b)  is  the  diapycnal  diffusive  flux  across  the 
isopycnal  b.  Formal  derivations  of  this  relation  appear 
in  Garrett  et  al.  (1995)  and  in  Marshall  et  al.  (1999). 

D{b)  consists  of  many  mixing  processes  that  can  re- 
sult in  diffusive  flux  across  isopycnal  b.  There  will  be 
contributions  to  D from  diffusive  mixing  in  the  thermo- 
cline,  while  horizontal  diffusive  fluxes  due  to  mesoscale 
stirring  would  arise  near  the  ocean  surface,  and  tidal 
mixing  near  rough  bathymetry  for  the  isopycnals  that 
intersect  with  the  ocean  bottom  in  such  regions.  In 
so  far  as  the  balance  in  the  upper  ocean  is  concerned, 
if  the  control  volume  is  chosen  to  be  bounded  by  two 
isopycnals  b and  b + db  and  a control  surface  at  the  bot- 
tom of  the  mixed  layer  (Figure  4),  Garrett  et  al.  (1995) 
showed  that  for  negligible  horizontal  mixing  and  with 
vertical  isopycnals,  vertical  mixing  D(b)  makes  no  con- 
tribution and  the  net  diapycnal  advection  is  only  due 
to  the  air-sea  transformation,  that  is, 

A(b)  = F(b)  (3) 

where  F(b)  = — Bo^f  is  the  air-sea  transformation, 
with  the  negative  sign  due  to  the  convention  that  the 
diapycnal  advection  from  less  buoyant  to  more  buoyant 
water  mass  is  considered  positive  (reversing  the  signs  of 
Walin’s  convention  in  Figure  1).  This  is  an  important 
idealization,  since  in  this  case  the  diapycnal  advection 
can  be  fully  diagnosed  by  surface  buoyancy  flux  and 
surface  buoyancy  outcrop  distribution. 

Speer  (1997)  evaluated  F(b)  and  A(b)  using  a hy- 
drographic section  at  11  °S,  as  shown  in  the  schematic 
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Figure  5.  A schematic  showing  the  control  volume  for 
the  diagnostic  calculation  carried  out  by  Speer  (1997). 


Figure  6.  The  ideal  algorithm  can  no  longer  be  ex- 
pected to  be  valid  if  time-dependent  effects  are  taken 
into  account  with  entrainment  at  the  mixed  layer  base. 


Figure  5.  These  diagnostic  calculations  showed  that 
the  residual  between  A(b)  and  F(b)  is  significant,  and 
therefore  transformation  due  to  the  diffusive  flux  term 
dD(b)/db  must  be  considered.  While  this  diagnostic 
calculation  cannot  be  used  to  pinpoint  which  processes 
contribute  to  the  D term,  Speer  (1997)  showed  that  the 
calculations  are  consistent  with  diapycnal  mixing  in  the 
thermocline. 

The  formulation  considered  thus  far  has  assumed 
steady  state,  but  an  accumulation  within  the  isopyc- 
nals  can  be  interpreted  as  a net  inflow,  and  (2)  can  still 
be  used.  A steady  state  scenario  neglects  the  role  of  di- 
apycnal advection  due  to  time  dependent  deepening  and 
shallowing  of  the  mixed  layer,  at  both  seasonal  and  di- 
urnal timescales.  Conceptual  modifications  are  required 
to  relax  the  steady  state  assumption  for  the  water  mass 
transformation  balance  as  shown  in  Figure  6. 

The  next  section  considers  diapycnal  advection  due 
to  the  ubiquitous  process  of  mixed  layer  deepening  and 
entrainment.  Upwelling  at  the  base  of  the  mixed  layer 


(a)  (b) 


Figure  7.  A schematic  showing  isopycnal  distribution 
(left)  and  buoyancy  distribution  (right).  The  solid  lines 
are  before  the  mixing  event  and  the  dashed  lines  are 
after  the  mixing  event. 


and  subsequent  mixing  results  in  water  mass  transfor- 
mation which  is  conceptually  very  similar  to  that  due  to 
mixed  layer  deepening  and  entrainment  at  the  base  of 
the  mixed  layer.  Regions  of  significant  upwelling  such  as 
those  found  near  submesoscale  instabilities  near  fronts 
are  considered  next  via  an  example  from  FASINEX 
( Weller  1991).  The  final  section  considers  whether  the 
stirring  due  to  mesoscale  processes  should  be  considered 
as  an  independent  process  for  water  mass  formation  and 
hence  considered  to  be  independent  of  the  transforma- 
tion due  to  the  mixed  layer  entrainment,  or,  do  the  two 
processes  happen  together,  in  series  and  are  interdepen- 
dent. The  water  mass  transformation  rates  would  only 
be  additive  in  the  former  case.  We  conclude  with  a 
summary  and  unanswered  questions. 


2.  Transformation  Due  to  Mixed  Layer 
Deepening 


Garrett  and  Tandon  (1997)  show  that  buoyancy  re- 
distribution in  the  vertical  by  mixing  achieves  diapycnal 
advection  across  an  isopycnal  surface.  Figure  7 shows 
an  isopycnal  that  outcrops  at  position  xq  due  to  vertical 
mixing,  and  the  buoyancy  balance  is  shown  on  the  right. 
x is  taken  to  be  the  direction  normal  to  the  outcrop 
towards  less  dense  waters. 

With  a uniform  buoyancy  gradient  in  the  vertical  and 
in  the  horizontal  everywhere,  the  volume  between  the 
successive  isopycnals  does  not  change,  i.e,  there  is  no 
diapycnal  advection.  However,  a departure  from  these 
conditions  implies  a net  diapycnal  advection.  In  par- 
ticular, if  a pre-existing  mixed  layer  of  depth  —h  + dh 
deepens  to  — h,  as  shown  in  Figure  8,  Garrett  and  Tan- 
don (1997)  show  that  the  net  diapycnal  advection  Fml 
is  given  by 


Fml 


-bxxb~3(Ab)2 


(4) 
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(a)  Isopycnal  position 


(b)  Buoyancy  profiles  at  x = x0 


b0b, 


-h  + 8h 

-h  --- 


bQ  — b]  ~ — Ab 


8h 


Figure  8.  A schematic  showing  isopycnal  distribu- 
tion (left)  and  buoyancy  distribution  (right)  for  a pre- 
existing mixed  layer  in  presence  of  both  bxx  and  bxz. 
The  solid  lines  are  before  the  mixing  event  and  the 
dashed  lines  are  after  the  mixing  event. 


per  unit  outcrop  length,  the  subscripts  i and  s refer 
to  the  gradients  below  the  mixed  layer  (interior)  and 
surface  mixed  layer  respectively,  Ab  is  the  buoyancy 
change  at  the  mixed  layer  base  before  a deepening  event, 
h is  the  mixed  layer  depth,  w is  the  vertical  velocity  at 
the  mixed  layer  base,  bxx,  bx  are  buoyancy  derivatives 
calculated  at  the  surface,  t is  the  time  and  T = 1 yr. 
Fml  is  the  integrated  diapycnal  advection  over  the  year 
for  each  deepening  event  w+  — >0,  since  the  upwelling 
with  vertical  velocity  w has  an  effect  similar  to  mixing 
and  subsequent  deepening  by  w St. 

The  Fml  expression  above  consists  of  two  terms,  the 
first  depending  on  vertical  variation  of  horizontal  buoy- 
ancy gradients  and  the  second  depending  on  the  cur- 
vature or  the  horizontal  variation  of  the  surface  buoy- 
ancy gradients.  Garrett  and  Tandon  (1997)  remark  on 
several  cases  in  the  limit  as  bxi,  the  buoyancy  gradi- 
ent at  the  mixed  layer  interface  and  bxs,  the  buoy- 
ancy gradient  at  the  surface  vanish,  showing  that  the 
singularity  in  (4)  is  no  worse  than  the  singularity  in 
Fair-sea  = — Bo(dS/db).The  formulae  need  information 
on  both  large-scale  horizontal  gradients  and  vertical  cy- 
cling of  the  mixed  layer. 

Tandon  and  Zahariev  (2001)  have  used  the  Marine 
Light  Mixed  Layer  (MLML)  experiment  mooring  obser- 
vations ( Plueddemann  et  al.  1995)  which  include  both 
spring  and  fall  mixed  layer  transitions  and  used  a com- 
bination of  mixed  layer  model  results  and  hydrography 
( daSilva  et  al.  1994)  to  calculate  the  bxx  terms  to  get 
an  order  of  magnitude  estimate  of  Fml-  The  sensitivity 
to  averaging  of  synoptic  events  is  also  explored.  Their 
calculations  indicate  that  if  hourly  winds  are  used,  the 
water  mass  transformation  due  to  mixed  layer  entrain- 
ment has  annual  peak  contributions  of  about  0(4)  Sv 
for  ert  = 24.0.  This  is  comparable  to  the  annual  trans- 
formation attained  by  diapycnal  mixing  in  the  upper 
ocean  water  masses  by  Zhang  and  Talley  (1998).  How- 
ever, with  daily  averaged  winds  and  without  diurnal 


Figure  9.  With  meanders  superposed  on  a mean  isopy- 
cnal distribution,  Fml  estimates  are  higher,  both  due 
to  the  variation  in  buoyancy  gradient  and  the  increased 
length  of  the  isopycnal  over  which  water  masses  get 
transformed. 


variation  in  buoyancy  forcing,  this  contribution  is  up  to 
an  order  of  magnitude  smaller.  Another  set  of  mixed 
layer  simulations  includes  an  annual  cycle  with  a shal- 
low and  strong  summer  thermocline.  Inclusion  of  syn- 
optic summer  forcing  for  this  scenario  leads  to  trans- 
formation values  several  times  larger  than  above,  about 
0(14)  Sv  at  ex t = 24.0.  The  peak  contribution  in  this 
case  is  almost  two  orders  of  magnitude  smaller  if  the 
synoptic  forcing  is  averaged  daily  and  the  diurnal  cycle 
is  not  resolved.  The  upper  bound  for  the  bxz  term  at 
the  MLML  experiment  site  is  0(3)Sv.  These  estimates 
are  likely  to  be  underestimates  since  hydrographic  data 
average  over  isopycnal  meanders,  due  to  both  increase 
in  the  curvature  terms  and  an  increase  in  the  isopycnal 
length  (Figure  9).  More  details  are  presented  in  the  ap- 
pendix of  Tandon  and  Zahariev  (2001).  Contribution 
due  to  mesoscale  and  sub-mesoscale  meanders  is  also 
conceptually  considered  later  in  more  detail  in  section 
4.2  of  this  paper. 

3.  Transformation  Due  to  Upwelling 
Near  Fronts 

While  the  bxz  term  was  found  to  be  small  in  Tandon 
and  Zahariev  (2001)  away  from  the  fronts,  this  need  not 
be  the  case  near  fronts.  Mesoscale  eddies  (wavelength 
~ 500  km)  strain  the  surface  into  0(100  km)  tongues, 
which  leads  to  sub-mesoscale  features  of  20-50  km  di- 
ameter which  induce  vertical  velocities  of  0(40  m/day) 
thereafter  (e.g.  Pollard  and  Regier  (1992)  for  Frontal 
Air-Sea  Interaction  Experiment  site,  Rudnick  (1996) 
and  Rudnick  and  Luyten  (1996)  for  the  Azores  front). 
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Due  to  the  rapid  upwelling  near  the  fronts,  the  up- 
welling  at  the  base  of  the  mixed  layer  fw  dt>Q  w dt  dom- 
inates the  deepening  Jdh>Qdh  in  the  Fml  term.  Using 
typical  FASINEX  values  of  bxs  = 1.5  x 10-7  s-2  (for 
temperature  gradient  of  IK/20  km),  bxx  = 1.5  x 10~12 
m-1s-2  and  A b = 10-3  m s-2.  With  upwelling  of  0(40 
m/day)  along  half  the  length  of  a convoluted  frontal 
length  of  3000  km,  the  annualized  transformation  val- 
ues are  0(0.2  Sv)  for  the  bxx  term  and  0(5  Sv)  for  the 
bxz  term.  The  Fml  term  is  therefore  likely  to  have  a 
greater  contribution  from  bxz  term  and  upwelling  near 
sub-mesoscale  features. 

To  summarize,  while  Fml  estimates  are  0(5)Sv,  they 
are  particularly  sensitive  to  the  resolution  of  synoptic 
events  and  sub-surface  density  structure  near  fronts. 
It  is  also  reasonable  to  ask  whether  these  effects  are 
already  assumed  to  be  present  in  a horizontal  diffu- 
sive flux  due  to  mesoscale  stirring.  Inspired  by  Gar- 
rett’s review  (this  volume)  as  well  as  Nakamura  (1996) 
and  Winters  and  D’Asaro  (1996),  we  next  ask,  under 
what  circumstance  can  these  two  processes  result  in  wa- 
ter mass  transformation  independent  of  each  other  and 
when  must  they  be  considered  interdependent. 

4.  Mesoscale  Stirring  and 
Mesoscale-Mixed  Layer  Interaction 

Horizontal  stirring  due  to  mesoscale  eddies  and  sub- 
sequent mixing  would  contribute  to  D(b),  and  sub- 
mesoscale  upwelling  contributes  to  Fm L(b).  How  should 
the  distinction  be  made  between  water  mass  transfor- 
mation due  to  these  two  processes?  How  do  we  ensure 
that  the  transformation  isn’t  being  diagnosed  multiple 
times  under  different  processes  which  are  in  series  (and 
hence  not  independent?  These  questions  are  clearly  re- 
lated to  the  how  the  resolved  and  unresolved  processes 
are  defined. 

Recent  understanding  ( Speer  et  al.  1997;  Garrett 
and  Tandon  1997;  Marshall  et  al.  1999,  Nurser  et 
al.  1999)  has  concentrated  on  expanding  (2)  heuristi- 
cally  in  terms  of  the  physical  processes,  which  can  be 
written  as 

A{b)  = Ask  man  (b)  + Aeddy(fr)  + ^mean  (b)  (5) 

rp  tp  dDe  ddy  d-Dint 

— ^ air —sea  < ^ ML  ^ ^ (o) 

with  (5)  for  the  dynamics  which  must  balance  (6)  for 
the  thermodynamics.  The  dynamical  processes  are 
the  diapycnal  volume  flux  due  to  wind  forcing  (Ek- 
man),  mesoscale  diapycnal  volume  flux,  and  the  mean 
geostrophic  flow  across  isopycnals.  The  thermodynamic 
processes  are  the  air-sea  transformation  Fair_sea,  the 
contribution  due  to  mixed  layer  entrainment  Fml,  hor- 
izontal mixing  due  to  eddies  dD^y/db,  and  the  interior 


z 


Figure  10.  An  illustration  showing  the  modified  con- 
trol volume.  A.Ekman  (b)  refers  to  the  wind-driven  diapy- 
cnal volume  flux  past  the  control  volume  surface  b , while 
AboIus  (b)  or  Aeddy  (5)  refers  to  the  diapycnal  volume  flux 
due  to  the  mesoscale  eddies  past  this  control  surface. 
The  lower  control  surface  C is  the  base  of  the  winter 
mixed  layer.  Dashed  vertical  fines  show  the  lateral  con- 
trol surfaces  for  the  dynamic  method;  the  instantaneous 
isopycnal  positions  shown  as  solid  fines  form  the  lateral 
control  surfaces  for  the  thermodynamic  method. 

diapycnal  mixing  dDint/db.  This  procedure  is  naturally 
isopycnic/diapycnic  and  allows  comparison  of  these  di- 
verse physical  processes  under  a single  framework. 

The  terms  that  are  most  readily  calculated  from  cli- 
matological data  in  the  above  framework  are  Aekman 
and  Fair_Sea,  e.g.  the  recent  evaluations  of  AEkman  and 
Fair_Sea  for  the  Southern  Ocean  by  Speer  et  al.  (1997). 
The  evaluation  of  Aekman  is  generally  done  follow- 
ing isopycnal  surfaces  around  the  annual  cycle.  How- 
ever, the  control  volumes  for  the  dynamic  and  thermo- 
dynamic methods  can  be  chosen  differently  as  long  as 
they  coincide  at  the  winter  mixed  layer  depth.  Thus, 
a modified  control  volume  for  the  dynamic  calculation 
can  be  chosen  ( Garrett  and  Tandon  1997)  such  that  this 
consists  of  the  volume  between  the  vertical  projections 
of  the  isopycnal  surfaces  below  the  winter  mixed  layer 
to  the  surface  (Figure  10).  Therefore,  this  modified  con- 
trolled volume  does  not  change  much  during  the  annual 
cycle. 

4.1.  The  triple  decomposition  into  mean,  eddy, 
and  turbulent  flow 

While  choosing  a fixed  control  volume  as  described 
above  simplifies  the  diagnostic  evaluations  of  the  dy- 
namic terms  in  (5),  the  issues  concerning  the  decompo- 
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sition  and  overlap  of  the  thermodynamic  terms  in  (6) 
still  remain.  To  clarify  the  role  of  horizontal  stirring 
dDn/db  and  mixed  layer  entrainment  Fml  , following 
Davis  (1994)  and  Garrett  (this  volume),  we  apply  the 
triple  decomposition  to  the  water  mass  transforma- 
tion ideas.  As  these  authors  stress,  such  a decompo- 
sition necessarily  presupposes  a gap  in  time  (or  space) 
scales  such  that  a three-way  decomposition  can  be  de- 
fined. For  the  water  mass  transformation  diagnostics 
the  buoyancy  field  is  split  into  the  mean,  eddy,  and  tur- 
bulent parts.  The  diapycnal  advection  un  is  also  split 
into  three  parts,  each  representing  the  mean,  eddy,  and 
turbulent  advection  (subscripts  m,e,  and  t respectively) 
across  the  mean  buoyancy  surface,  i.e, 

b = bm  + be  + bt  and  (7) 

Wn(^m)  = ’U'nm^Prri)  ~b  ^ti  ei^m  ) + unt(b  m )•  (8) 

Here  we  envisage  the  mean  buoyancy  field  bm  to  be 
the  modified  mean  buoyancy  field  ( Killworth  2001,  Mc- 
Dougall  and  McIntosh  2001)  although  for  the  large  scale 
water  mass  transformation  diagnostics  the  difference 
with  Eulerian  mean  buoyancy  field  is  probably  small. 
We  denote  ( ) for  averaging  over  a time  scale  large  com- 
pared to  turbulence  but  small  compared  to  that  of  the 
eddies  and  < ( ) > for  averaging  over  a time  scale  large 
compared  to  that  of  the  eddies  but  short  compared  to 
that  of  the  mean.  Substituting  into  the  instantaneous 
buoyancy  equation  and  subsequent  averaging  over  the 
turbulent  time-scales  and  the  eddy  time  scale  yields 

dbm  dbe  dbt 

Unm ' "a^+  < Une ' >_  ~ < Unt ' > (9) 

if  the  mean  flow  is  considered  steady  at  long  time-scales. 
The  right  hand  side  of  eq.(9)  is  the  — term  and  the 
terms  on  the  left  hand  side  are  mean  and  eddy  advec- 
tion of  mean  and  eddying  isopycnals  respectively.  This 
equation  is  subject  to  the  surface  boundary  condition 

< webe  > + < wtbt  >=  Bo  (10) 

thus  suggesting  that  the  surface  buoyancy  flux  is  bal- 
anced by  both  the  vertical  buoyancy  flux  due  to  eddies 
and  turbulent  flux  at  the  surface.  Garrett  (this  volume) 
discusses  the  buoyancy  variance  budget  for  this  decom- 
position. 

4.2.  Mesoscale  transformation  and  the 
entrainment  contribution:  In  series  or  parallel? 

If  the  diagnostics  are  done  in  a framework  that  does 
not  resolve  eddies,  then  the  transformation  relationship 
can  be  written  in  terms  of  buoyancy  balance  about  bm  : 

A(6m)  = Amean(^m)  “b  Aeddy(^m)  "b  Ae]c  man  (^m  )(n) 


(12) 

where  the  contribution  ^*(6m)  is  due  to  the  conver- 
gence of  the  eddy  stirring  horizontal  diffusive  flux  and 
FML{bm)  is  the  contribution  to  water  mass  transfor- 
mation due  to  mixed  layer  deepening  and  entrainment, 
both  evaluated  about  the  mean  isopycnal  positions. 

Following  the  triple  decomposition  suggested  above, 
if  the  eddies  are  resolved  in  this  framework,  the  water 
mass  transformation  balance  becomes 

A(bm^  — Amean(6m)  -{-  Aeddy(^m)  “b  AeIc  man(^m)  (13) 

= F(bm  + be)-FuL(bm  + be)-~(bm  + be) 

(14) 

and  the  eddy-stirring  component  does  not  arise. 
The  contributions  in  eq.  (11)  and  eq.  (13)  are  the  same 
due  to  the  choice  of  control  volume  in  Figure  10.  How- 
ever, the  contributions  in  eq.  (12)  and  eq.  (14)  are 
different.  In  this  case,  the  diapycnal  advection  due  to 
horizontal  mixing  by  the  eddies  and  the  contribution 
due  to  mixed  layer  entrainment  are  in  parallel  if  the 
budgets  are  considered  about  the  mean  buoyancy  sur- 
faces ( bm ) that  do  not  resolve  the  eddies  as  in  eq.  (12). 
However,  these  two  processes  are  in  series  if  budgets 
are  considered  for  eddying  buoyancy  surfaces  ( bm  + be) 
as  in  eq.  (14),  where  eddying  isopycnals  are  defined  by 
modified  mean  density  that  averages  over  a time  scale 
long  compared  to  the  turbulence  but  short  compared  to 
the  eddies.  Thus  the  Fml  estimates  computed  by  Tan- 
don  and  Zahariev  (2001)  (section  2)  which  are  based  on 
mean  monthly  isopycnals  that  average  over  the  eddies 
correspond  to  the  estimates  for  FW/,(6m),  while  the  es- 
timates earlier  near  frontal  zones  (section  3)  correspond 
to  the  latter  case  FML(bm  + be). 

5.  Discussion 

The  uncertainty  in  Fml  remains  unacceptably  high 
in  both  model  based  estimates  of  coarse  non-eddy  re- 
solving simulations  (e.g.  Nurser  et  al.  1999)  and  data 
based  evaluations  which  are  based  on  localized  mixed 
layer  simulations  ( Tandon  and  Zahariev  2001).  More 
accurate  basin  wide  estimates  are  necessary  to  make 
progress.  While  considering  budgets  in  an  eddy  resolv- 
ing framework,  the  water  mass  transformation  due  to 
eddies  happens  in  series  with  diapycnal  mixing  and  en- 
trainment. In  this  case,  another  significant  concern  is 
identifying  the  rate  controlling  process  for  diapycnal 
transformation  in  the  upper  ocean — do  the  mesoscale 
eddies  define  the  rate  at  which  net  diapycnal  advection 
will  happen  or  does  it  primarily  depend  on  the  rate  at 
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which  diapycnal  mixing  occurs?  The  mesoscale  stirring 
would  be  dominant  for  water  mass  transformation  if 
the  mixing  processes  are  so  efficient  in  the  upper  ocean 
that  they  act  to  dissipate  the  eddies  quickly  and  trans- 
form water  masses  on  a fast  time  scale.  In  this  case, 
it  is  the  baroclinic  instability  and  larger  scales  that 
would  set  the  net  diapycnal  advection  rates.  On  the 
other  hand  it  is  possible  that  the  mixed  layer  processes 
allow  non-linear  interactions  amongst  eddies,  and  the 
net  diapycnal  advection  is  then  significantly  dependent 
on  mixing  processes  in  the  upper  ocean.  While  recent 
eddy  resolving  simulations  have  concentrated  on  adi- 
abatic properties  of  mesoscale  paramaterizations  (e.g., 
articles  by  Marshall,  Killworth,  and  McDougall  from 
this  volume)  diabatic  interactions  in  the  upper  ocean  re- 
main unexplored,  and  the  question  above  remains  unan- 
swered. Determining  which  process  is  dominant  for  wa- 
ter mass  transformation  rate  should  be  feasible  with 
multiple  eddy-resolving  simulations  whose  sensitivity  to 
mixing  in  the  upper  ocean  and  air-sea  interaction  is  ex- 
plored systematically. 
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